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g : ABSTRACT 

o 

^ ' We study the FCNC top quark decays t ^ ch in the framework of the MSSM, 

where h = h^, H^, A^ is any of the supersymmetric neutral Higgs bosons. We 
include the leading set of SUSY-QCD and SUSY electroweak contributions. 
, While the FCNC top quark decay into the SM Higgs boson has such a negli- 

(D I gible rate that will not be accessible to any presently conceivable accelerator, 

1-^ ' we find that there is a chance that the potential rates in the MSSM can be 

measured at the high luminosity colliders round the corner, especially at the 

^ . LHC and possibly at a future LC, but we deem it difficult at the upgraded 

H I Tevatron. In view of the large SUSY-QCD effects that we find in the Higgs 

channels, and due to some discrepancies in the literature, we have revisited 
the FCNC top quark decay into gluon, t ^ eg, in our framework. We con- 
firm that the possibility of sizeable rates does not necessarily require a general 
pattern of gluino-mediated FCNC interactions affecting both the LH and the 
RH sfermion sectors - the LH one being sufficient. However, given the present 
bounds on sparticle masses, the gluon channel turns out to lie just below the 
expected experimental sensibility, so our general conclusion is that the Higgs 
channels t —>■ ch (especially the one for the light CP-even Higgs) have the 
largest potential top quark FCNC rates in the MSSM, namely of order 10~^. 
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1 Introduction 

The study of the virtual effects in top quark decays into Higgs bosons could be the clue 
to physics beyond the Standard Model (SM). This fact has already been demonstrated 
for top quark decays into charged Higgs bosons, t -^ H^h [Q , within the context of the 
Minimal Supersymmetric Standard Model (MSSM)p|0. Indeed, the potential existence of 
large quantum effects induced by supersymmetric particles in certain regions of parameter 
space may lead to highly significant changes in the partial width of that decay, and this 
feature could have a serious impact on the Higgs searches at the Tevatron 0, |[ • A 
situation which is in contrast to the SM decay of the top quark, t -^ W^ b, where the 
SUSY quantum effects are in general much more modest 0. 

Similarly, one may expect that the top quark decays into the neutral Higgs bosons of 
the MSSM may undergo relevant enhancements. Notice, however, that in contradistinc- 
tion to the charged current decays mentioned above, the loop contributions are in this case 
the lowest order effects as the neutral Higgs decays of the top quark are mediated by Fla- 
vor Changing Neutral Currents (FCNC). Therefore, since the FCNC processes are rather 
suppressed the possibility of MSSM enhancements should be very welcome, especially for 
the physics program at the LHC and perhaps also at a future linear collider LC. In the 
LHC, for example, the production of top quark pairs will be very high: a(tt) = 800 pb 
- roughly two orders of magnitude larger than that of the Tevatron II at i/i = 2 TeV. 
In the so-called low-luminosity phase (10^^ cm~^s~^) of the LHC one expects about one 
tt-pair per second, that is to say of the order of ten million tt-pairs per year [0]. And 
this number will be augmented by one order of magnitude in the high-luminosity phase 
(lO^^cr/i^^s^^). As for a future e"*"e~ linear collider running at e.g. ^/s = 500 GeV, 
one has a smaller cross-section a{tt) = 650 fb but a higher luminosity factor ranging 
from 5 X 10'^^ cm~'^s~^ to 5 x 10^^ cm"'^s^^ [§] and of course a much cleaner environment. 
One thus expects that both the LHC and the LC will initially deliver datasets of order 
10/6~^/year increasing to several 100/6~^/year in the high-luminosity phase. On the 
other hand, at the Tevatron II during the highest luminosity era (TeV33) one expects 
typical datasets of 30/6"^. It follows that if the branching ratios of the FCNC decays are 
augmented by extra contributions beyond the SM, one should be able to collect enough 
statistics (perhaps some few hundred to few thousand events) from the combined output 
of these machines enabling us to perform an efficient study of these rare decays. 

We should immediately point out that, within the strict context of the SM, the pos- 
sibility of detecting FCNC decays of the top quark is essentially hopeless. In particular, 
it has recently been recognized that the FCNC rate of the top quark decay into the SM 
Higgs boson (t -^ cHsm ) is much smaller p, |^ than originally thought |ll[]: It turns 
out that BR{t -^ cHsm) = 1 ■ 10"^^ - 4 ■ lO'^^ {Mz < Mh < 2 Mw) §, which means 
that it is far out of the range to be covered by any presently conceivable high luminosity 
machine. On the other hand, the situation with the FCNC decays of the top quark into 
gauge bosons (t — >■ c\^; V^ = 7, Z, g) is not much more promising in the SM, the branching 
ratios being at most of order 10~^^ for the photon, slightly above 10~^^ for the Z-boson, 
and at most 10~^° for the gluon channel Jll], |12[. Thus the highest SM rate, namely 
that oi t ^f eg, is still 5 (resp. 7) orders of magnitude below the feasible experimental 
possibihties at the LHC (resp. Tevatron II). Clearly, detection of FCNC decays of the 
top quark at visible levels (viz. BR{t — * cX) > 10~^ — 10~^) by any of the future high 

^The corresponding study for general two-Higgs-doublet models (2HDM's) is also available in Ref. |3|. 



luminosity colliders would be instant evidence of new physics! 

Therefore, one may judiciously ask whether extra virtual effects beyond the SM can 
help to bring the top quark FCNC decay ratios to within observable levels. For exam- 
ple, for the t ^ cV decays one finds, within the non-supersymmetric 2HDM's, that 
there can be significant enhancement s |TT], |T3| which, however, turn out to be insufficient. 



A similar situation occurs in the more interesting case of the MSSM, where again in 
spite of the potential enhancements the electroweak gauge boson channels fall short to 
be detected [0]-|jl3. Only the FCNC top quark decay in the gluon channel could be 
fairly sensitive to the SUSY corrections in non-negligible regions of the MSSM parameter 
space |J4[-|T^, a fact which we wish to revisit within our framework in order to compare 



with the more exceptional possibilities offered by the Higgs channels - on which we will 
mainly concentrate^. 

These Higgs channels comprise the FCNC top quark decays into the two CP-even 
( "scalar" ) states and the CP-odd ( "pseudoscalar" ) state of the Higgs sector of the MSSM 

0, 

t^ch {h = h'',H^,A''). (1) 

Worth emphasizing is the fact that in the MSSM (in contrast to the SM or the uncon- 
strained 2HDM) at least one of these decays (viz. t -^ ch^) is always possible, for in 
the MSSM there is an upper bound on the mass of the lightest CP-even Higgs boson, 
M/iO < 135Gel^||2y]-|2^, which is below the top quark mass. Moreover, for a sufficiently 



light pseudoscalar mass, M^o < mj, all three decays (^ are in principle possible, if the 
SUSY masses are not that high so as to induce too large positive corrections to M^^. 

As it is also the case with their charged- current counterpart mentioned in the begin- 
ning, these decays could be greatly enhanced in wide regions of the parameter space. 
Although some work already exists in the literature on FCNC decays of the top quark 



into Higgs bosons within the framework of the MSSM |2J] , we feel that it is still rather 
incomplete because it fails to include some of the most significant contributions and it 
does not make use of the one-loop Higgs mass relations of the MSSM [|^-[^. These 



relations play an essential role in correlating the various quantum effects for the different 
channels (P and are fundamental in establishing the aforementioned upper bound on 
the mass of the lightest supersymmetric Higgs particle. Therefore, our purpose is to go 
beyond those preliminary calculations and show, from a more rigorous and systematic 
treatment of the different kinds of quantum effects and of the MSSM parameter space, 
that the FCNC width of the top quark could in fact reach the experimentally visible level 
at the high-luminosity colliders |^ . 



The paper is organized as follows. In Sect. 2 we give the SUSY Lagrangian interactions 
relevant for the FCNC decays of the top quark. In Sect. 3 we report on the SUSY 
electroweak (SUSY-EW) and supersymmetric QCD (SUSY-QCD) one-loop contributions 
to the t ^ c h decays (|l]). In Sect. 4 we address the SUSY-QCD contributions to t ^ c g 
in our framework and compare with t ^ c h. Finally, in Sect. 5 we further discuss our 
results and deliver our conclusions. 



^For recent studies on FCNC top quark decays beyond the MSSM, see Ref. |18| and references therein. 



2 Relevant interaction Lagrangians 

The MSSM interaction Lagrangian involving fermions and SU{3)c x SU{2)l x U{1)y 
gauge bosons is well-known, and will not be spelled out here in any detail 0, |^. We will 
just focus on the relevant interaction pieces entering the main one-loop contributions to 
the Higgs channels (|l]). 

To better present the analytic results of this computation we shall define first a nota- 
tion that allows us to treat the three possible decays (|l]) in an unified way. We introduce 
a vector array of neutral Higgs fields 

$° = (if °, /i°, A°) , (2) 

and another one for the charged Higgs and Goldstone bosons 

<!?+ = {H+,G+) . (3) 

Then the interaction Lagrangian of quarks with neutral and charged Higgs bosons in the 
MSSM can be written 



grriu 



^*.. = - oJ 1^ E <^'MK'r''PL + {Kn'PR) 



2Mwsmp ^±^ 



u 
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gmd ^^-^ /f.OJ/T^Odr> , /T^Od\*: 




+ ^^^"'^ E [^^ ^(Kl^Pl + K^'Pr) u + h.c] , (4) 

\/lMw r-=l,2 

with Pl,r = (1/2) (1 =F 75) the chirality projectors. We have defined the following sets of 
i^-matrices: for the neutral Higgs sector 



ir°"= cosa , K^.'^= I -sina | , (5) 

\ i cos (3 J 

and for the charged Higgs sector 

K^l' = m, ( ^^_^/ ) , K^r = m„ ( '°^* ^ ) . (6) 

In the above formula Vud is the CKM matrix element (assumed to be the same as for the 
W^ boson interactions with ordinary quarks). 

The relevant Yukawa couplings involving charginos, quarks and squarks are contained 
in 

^nJx+ = -9 Vud dl ^t (^Vt? Pl + A^':^ Pr) u + h.c. , (7) 

with u [d) up-type quarks (down-type squarks) of any generation, the coupling matrices 
being ^ 

where 

" ~ ^ ~ V2Mw sin/3 ' ^~ 9 ~ y^Mw cos/5 ' ^ ' 
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are the up-like and down-like Yukawa couplings normalized with respect to the SU{2)l 
gauge coupling. Of fundamental importance for the SUSY-EW enhancement of our FCNC 
decays through these Yukawa couplings is the value of the parameter tan/3 = f2/fi[|l9|. 
In the equations above R^'^^ is the 2x2 matrix that diagonalizes the squark mass squared 
matrix in chiral space through 

?1 = E^&' (10) 

b 

q'^ being the weak-eigenstates and qa the mass-eigenstates. 

Finally (regarding the EW part) we quote the interaction Lagrangian for triplet Higgs 
vertices in the MSSMQ 

U^^ = -9Y. Brst $+ $; $? , (11) 

r,s,t 

and the chargino couplings to neutral Higgs bosons 

^x+x+* = -9T.^'rXtm^LPL + Wl^^Pn)xi . (12) 

In both cases we have encapsulated the remaining notation in two 3x3 matrices Brst 



and Wlj whose explicit form can be identified from Ref. |T9|. These matrices just give the 
corresponding Feynman rules (divided my —ig). 

On the other hand, the necessary SUSY-QCD interactions (in the mass-eigenstate 
basis) for our FCNC decays are contained in the Lagrangian 



9_ 

V2 
9. 

V2 
9s 

V2 



-CSUSY-QCD = -f^^c [-RL Pl - Rla Pr] (la,i Kj tj 

^ Vf [Rla Pl - Rla Pr] C. a?,. 9 



^f [Rl^ Pl - i?L Pr] Ta, K, n, + h.c. , (13) 



where tp^ stands for the gluino spinor and A^ are the SU{2))c Gell-Mann matrices. More- 
over, the 6x6 rotation matrices i?^'^-* generalize those in eq. ([I0| ) and are needed to 
diagonalize the squark mass matrices in (flavor) x (chiral) space as follow^: 

R^'^^^MlR = Mlj, = diag{ml,...,ml} {q^u,d), (14) 

where tM?- t-, is the 6x6 square mass matrix for squarks in the EW basis {q'a), with indices 

a = 1, 2, 3, . . . , 6 = ul, uji, cl, . . . ,iji for up-type squarks, and a similar assignment for 
down-type squarks. However, in this study only up-type squarks are involved, so that we 
understand that the above diagonalizing matrices refer to them. The intergenerational 
mixing terms leading to gluino-mediated FCNC couplings lie in the off-diagonal entries 
of the mass matrices. However, in order to prevent the number of parameters from being 

■^Note that the elements Brss are complex and Bsss — 0. 

^First latin indices (a, b, . . .) refer to different squark states of the same flavor and take values 1, 2, with 
rriq-^ < rriq^ in the mass-eigenstate basis and q',^ ^x — 'i'il r} ^^ ^^^ weak-eigenstate basis. Latin indices 
i,j,... — 1,2 refer to charginos. r,s, . . . refer to Higgs and Goldstone particles. Greek indices a, (3, . . . 
refer to squarks in the 6x6 (flavor) x (chiral) space and take values 1, 2, . . . , 6, rriq-^ < rriq^ < . . . < rriq^ . 



too large, we have allowed (symmetric) mixing mass terms only for the left-handed (LH) 
squarks. This simplification is often used in the MSSM and it is justified by Renormal- 
ization Group (RG) analysis |]2^ . Following this practice, we introduce intergenerational 



(that is to say, flavor mixing) coefficients 6ij in the LL block of the mass matrix (namely 
the one involving only LH fields of any flavor) as follows: 

{MlL)ij = ml = 6ij rrii rrij {i ^ j) , (15) 

where rrii is the mass of the left-handed ith squark, and m?- is the mixing mass matrix 
element between generations i and j. Therefore, if the coefficients Sij are non- vanishing, 
for some i j^ j, then the structure of the diagonalizing matrices R^'^^ defined above must 
necessarily lead to gluino-mediated tree-level FCNC between quarks and squarks in the 
SUSY-QCD Lagrangian (p!3|). This scenario can be generalized if we further introduce 



FCNC interactions on the right-handed (RH) block of the mass matrix|2g, ^ (see Sects. 
4-5). 

We note that the induced FCNC couplings in the SUSY-QCD Lagrangian (p!3D ulti- 
mately stem from the fact that the squark mass matrix in general need not to diagonalize 
with the same matrices as the quark mass matrix, i.e. the so-called misalignment of quark 
and squark mass matrices. 

3 SUSY contributions to t ^ c h 

From the previous interaction Lagrangians, the computation of FCNC processes at one 
loop in a renormalizable theory is straightforward. In fact, renormalization of parameters 
and Green functions is not needed, unlike the usual flavor- conserving processes; one just 
computes the different diagrams that contribute to the process and the flnal result ob- 
tained after adding up all the amplitudes must be finite since no lowest order interaction 
could absorb the left over infinities. 

For every decay process ([l|) it is convenient to define an "effective" interaction vertex 

-^T = -tg u,{p) {Fl Pl + Fr Pn) Ut{k) , (16) 

where F^ and Fr are form factors associated to each chirality projector. They follow from 
explicit calculation of vertices and mixed self-energies. 

The one- loop graphs for the decays under study are depicted in Figs. |l], || and H. On the 
one hand the contribution from the vertex diagrams (Figs. ^ and |^) is obtained by direct 
calculation. On the other hand for the (mixed) self-energy diagrams (Cf. Figs. ^ and ^d) 
it is convenient to define (in analogy with the charged-current case0) the following 
structure: 

S,,(A;) = |^S^(P)P^ + ^S^(P)P^ + m, (Sz..(P) P^ + Sr,(P) P^) . (17) 

Here the mt factor multiplying the scalar part is arbitrary and it is inserted there only to 
preserve the same dimensionality for the different Sj. From the Lagrangian interactions 
introduced in Sect.Q we find that the effects of the mixed self-energy diagrams to the 
amplitude oit^c^^. take on the following general form in terms of the various S^: 



_ _ -igmt 1 r 



-^T^s. 



+ PR{Kfr [L^R]}ut{k), 



-igrric 



rrit 



2Mvysin j3 mf — ml 
+mc (S^(mj) + Si 



u,{p)[PlK' [mi(SiK) + S 



fel"^i. 



m. 



P^(K°^)* [L^R]}ut{k) 



where we have called Tg. the amplitude corresponding to attaching the external rth Higgs 
particle line to the zth quark leg (r = 1,2,3; i = c, t). From (0) the corresponding 
contribution to the form factors in eq. (pIBf) is transparent. 



3.1 SUSY-EW effects 

We start by reporting on the SUSY-EW effects, which by definition include the sparticle 
contributions plus the Higgs-Goldstone boson diagrams computed in the Feynman gauge 
(Cf. Figs.[l| and |]). For the electroweak part we work in the so called Super-CKM basis, 
that is, we take the simplification that the squark mass matrix diagonalizes as the quark 
mass matrix, so that FCNC processes appear at one-loop only through the charged cur- 
rent interactions (from charged Higgs and charginos) and with the same mixing matrix 
elements as in the Standard Model (the CKM matrix). It means that for the electroweak 
effects we do not take into account the possible mismatch between the matrices diagonal- 
izing the squark sector and those diagonalizing the quark sector as this would only result 
in a subleading additional correction. 

Let us consider the computation of the contributions from diagrams in Figs. |1| and 0. 
The graphs in which FCNC are mediated by gauge bosons (IV''') have not been included 
since they are subdominant as compared to the leading effects from the Yukawa couplings 
of the top and (at large tan j3) bottom quarks. Still, the leading terms from longitudinal 
W^ are taken into account through the inclusion of Goldstone bosons in the Feynman 
gauge. 

The simplest SUSY-EW contributions are the self-energy diagrams (Fig. |]). They are 
obtained after substituting the following expressions into the generic eq. (|13): 



s«(p; 
Sl(^^; 

mt T.Rs{k'^] 
nit ^Ls{k'^] 



^Rik'] 



^dk' 



mt Sfe(/c^ 



(a) 
(a) 
(a) 
(a) 

(b) 
ib) 
ib) 



I g" Vu Kd ^^ A^''2 Biik, M„ mj 
I g" Vtd Vai Af:} Af2 BiiK M„ mj 



I / M, Vt, K, J^ ALT/ Bo{k, M,, m^^, 



ig^MiVtdVcdA 



{d,c) .{d,t) 
+ai —ai 

{d,c) .{d,t) 



Bo{k,Mi,m_ 



da' 



I g m^ mt 



2M^ 



VtdVr 



td Vcd 



+ {Bo 



2M^ 



Vt,Vr 



td Vcd 



cot'/? (5o + 5i) (A;, MH±,mrf) 
Bi){k,Mw,md) , 
tan'/? {Bo + B,){k,MH±,md) 



+ iBo + B,){k,Mw,md) 



i g"^ iTLt m\ 
2M^ 



Vtd Vcd [Bo{k, MH±,md) - Bo{k, Mw, rrid)] 



nit ^Ls{k^ 



(b) 



i g^ rric m'^ 
2M2, 



Vtd Vcd [Bo{k, MH±,md) - Bo{k, Mw, rud) 



(19) 



where the two-point functions Bi are as in Ref.[^. Of course in the above expressions 
we understand a sum over the free indices a,i = 1,2 and over the three down-quark 
generations {d = d,s,b). 

The contribution to the form factors Fl and Fr in eq. (|TB|) from SUSY-EW vertex 
diagrams is much more cumbersome. Diagrams (a) and (d) of Fig. |l] give a generic con- 
tribution of the form 



Fl 

Fr 



Na 



l(l) 4(2) 



(1) Am 



l(l) /I (2) 



(Ci2 - Cn)m, A'^' A'^' - C,^ m, A'l> A'l> + Co m^ A'^' A'i 



Fl{A 



(*) 



A 



(*)^ 



(20) 



whereas diagrams (b) and (c) have the general structure 



A^, 



D 



i(i) n(2) n(3) 



(1) n(2) n(3) 



Co (Di^^ ^r^ Z^}?^ m, mi + Di^^ D}:^^ Dt' m, moi 

+Z}iJ' d[') Df nit mm + /^ij^ /^f /^f m^i m^s) 
+Ci2 m, (Dg) Dg) Df m, + D« Z^f Dg) m, 

+D«DfDfm^, + D«DS?)Dfm^,) 
+ (Cn - C12) m, (Dg) Dg) Z^g) m, + D^^ ^g^ Z^g^ m, 

+Z^g) Z^g) Z^g) mn. + Z)g) D^^^ Z^g) m..) + Co Z^g^ Z^g^ D^^^ 



'R 



Fl{D 



(*) 



D 



(*)^ 



(21) 



where the three-point functions Cj and Cij are as in Ref.[|l], ^f. 

Specifically, each vertex diagram of Fig. ^ contributes the following to the process 



Diagram (a); make the following substitutions in eq. (120): 



(1) 



A 



niA 



1 {d,c) 



is;' 



-bi ' 



If 



4^*) 4(2) 



A 



id,t) 



M, , Na = I g^VtdV^R^:';! {RjiT G^ 



where G^g is the well-known Feynman rule |T9[ for the vertex $[!(i'r d'* divided by 
—ig, in the electroweak-eigenstate basisQ. In this diagram the various three-point 
functions on ( PDD must be evaluated with arguments 



''da ' df, ) 



C* = C^{k, — p', Mj, m 
The convention for the momenta can be seen in Fig. |l]a. 
Diagram (b); make the following substitutions in eq. (^ 



D 
D 



(1) 

L 

(3) 



D 



Ad,c) 
4(d,t) Y) 



(1) 
R. 

(3) 



A 
A 



{d,c) 

-a] 

{d,t) 



D 



(2) 



w. 



ijL ' 



L) 



(2) 
R 



w. 



ijR 1 



mz)i = Mi , mD2 = ^i , Nd = i g^ Vtd Vcd 
C, = C,ik,-p',ni^^,Mi,Mj), 



^We recall that our sign convention ||l|, g, |2^ for the /i parameter is opposite to that of |19) 



Diagram (c); substitute in eq. (PI|): 



Df = Kt^' , /^g^ = Kt^' , 



g^rud 



rriDi = mD2 = rrih , Nd= i-rj^ ^V^tdK^ , 

4M{y cos/5 

C* = C*(A;, -p, m^+,md, ma) , 
• Diagram (d); substitute in eq. (pOl) : 

/l(l) _ IV'+crf 4(1) _ TV' + cd 4(2) _ Ty'+td 4(2) _ ^+td 

mA = md , Na = T'TrTTT ^sw Vtd V^d , 
C^ = C^{k,-p',md,m^+,m^+) . 

As can be noted from the above expressions, the form factors induced by Higgs me- 
diated diagrams - see Figs. 0c and d - have the property Fl = Fr for H^ and /i°, and 
Fl = -Fr for A\ 

We have performed the usual checks of the computation, in particular we find that 
the form factors Fl and Fr are free of divergences before adding up the three quark 
generations, both analytically and numerically in the implementation of the code. 

3.2 SUSY-QCD effects 

Using the Lagrangians (H) and (0) one can find the SUSY-QCD one-loop contributions 
to the processes under study. They are much more simple than the electroweak ones. 
The Feynman diagrams are depicted in Figs.|^a and b. The one-loop mixed self-energy in 
Fig. 0b is determined from 

Sl(^^) = -i2nasCFR3aRtaBii-k,mg,muJ , 

Si?(A;^) = -i2 7iasCFR4aRlaBi{-k,mg,muJ , 

mtT.Lsik'^) = -i2'KasCFTngRiaRlaBo{-k,Tng,muJ , 

mtT.Rs{k^) = -i2 7iasCFmgR3aRlaBQ{-k,mg,muJ , (22) 

where Cf = {N^ — l)/2Nc = 4/3 is the quadratic Casimir invariant of the fundamental 
representation of SU{3)c- 

Finally, the SUSY-QCD vertex contributions to the form factors (|16]) follow from 
Fig. 0a, and read 

Fl = N [mt Rap -Rga (C*!! — C'12) + "^c -R3/3 -RL C'12 + ^ Rap -RL^'o] , 

Fr = Fz. (3 ^ 4 , 5 ^ 6) , 

N = i%TT a^C F R*^pGl^^ Rsa , 

C* = C^{-k,p',mg,mu^,mu^) , (23) 

where G^^ is the well-known Feynman rule (divided by —ig) |jl9| for the vertex $°{t' u'/ 
in the electroweak-eigenstate basis. From these expressions it should be clear that if the 
up-type squarks would be degenerate, then by the unitarity of the i?-matrices the gluino- 
mediated FCNC effects would vanish (GIM mechanism). We have used this analytical 
property as an additional check of our numerical code. 

9 



4 Numerical Analysis oi t ^ c h 



After squaring the matrix element (|T^), and multiplying by the phase space factor, one 
obtains the decay width oi t —* c h, 



T{t-^ch) 



with A(a;^,y^,z^) = 
to define the ratio 



9' 



32 7nnt 



\\l2i 2 2 2\ 



{m^^+ml-ml)i\FL\' + \FR\^) + 2mtm,iFLF*^ + FlFR) , (24) 



{x'^ — {y+zY){x'^—{y—zy) the usual Kallen function. It is conventional 



B{t^ch) 



T{t^ch) 



(25) 



ch) of the decay mode, as 



V{t-^hW+) 

which will be the main object of our numerical study. 

This ratio is not the total branching fraction BR{t 
there are many other channels that should be added up to the denominator of (^) in the 
MSSM, if kinematically allowed, such as the two and three body decays of the top quark 
into SUSY particles, and also the charged Higgs decay channel t -^ H^h [|T], |2U[. For the 
mass spectrum used in our numerical analysis the former decays are phase space closed, 
whereas the latter could have a sizeable branching ratio. However, for better comparison 
with previous analyses of FCNC top quark decays pf-ljl^, the fiducial quantity ( p5D should 
suffice to assess the experimental viability of the FCNC decays under consideration. 

In the following we will analyze the numerical contributions to (^) from the SUSY- 
EW and SUSY-QCD sectors, in a sparticle mass model motivated by model building and 
RG analysis. However, we do not restrict ourselves to the spectrum of specific SUSY- 
GUT models - such as SUGRA or constrained MSSM modelsp], |3^. Furthermore, as 
announced in the beginning, all over our numerical analysis we use the full structure of 



the one-loop relations for the parameters in the Higgs sector of the MSSM |21 



We start with the EW effects. Although we have generally scanned the MSSM param- 
eter space up to 1 TeV level, the following input set has been chosen where the individual 
parameters have to be fixed at particular values to illustrate our results (as in Fig. H): 



tan /3 = 35 , /i = 


-500^61^ , 


M = 150 GeV , 


Mao = 100 GeV , 


mi^ = 150 GeV , 


m^^ = rrig = 


-- 200 GeV , At = 


-- Ag = 300 GeV , Ab 



-300 GeV . 

(26) 
We have taken the third generation quark masses as rrit = 175 GeV and rrib = 5 GeV. 
In the previous equation mi^ , mi are the lightest t and b mass, and all the masses are 
above present experimental bounds PT|. However, we have chosen a SUSY mass spectrum 
around 200 GeV, which is not too light, so the results will not be artificially optimized. We 
have also checked that all through our numerical analysis other bounds on experimental 
parameters (such as 6p) are fulfilled. 

We have found that the contributions to the form factors of (^) are of the same 
order for the chargino (Figs.|l|a,b and ^) and Higgs particles ||3^ (Figs.|I]c,d and ^) 

p^ . It turns out that they can be either of the 
The chosen negative value for Ab is to make the two 



-not included in previous analyses 
same sign, or of opposite sign. 



contributions of the same sign. In both cases Fr ^ Fl- This can be easily understood 
by looking at the interaction vertices involving (charged) Higgs-bottom-charm |T^ and 
chargino-sbottom-charm, where the latter can be tracked down from the explicit structure 
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of the Lagrangian (|^). In both of them the contribution to the right-handed form factor 
will be enhanced by the Yukawa coupling of the bottom quark, whereas the charm Yukawa 
coupling contributes to the left-handed form factor. We have checked that the inclusion 
of the first two generations of quarks and squarks only has an effect of a few percent on 
the total result. 

In Fig.^ we can see the evolution of the ratio (pSj) with various parameters of the 
MSSM by taking into account only the electroweak contributions. The growing of the 
width with tan/3 (Cf. Fig.^a) makes evident that the bottom quark Yukawa coupling 
plays a central role in these contributions. The evolution with the trilinear coupling At, 
-the main parameter that appears in the interaction vertex b^ &_r h- displayed in Fig. ^ 
shows that this parameter can enhance the width by some orders of magnitude. On the 
other hand the evolution with he higgsino mass parameter fi (Cf. Fig. ^) is comparatively 
mild in the region away from the origin /i = 0. The shaded region in Fig. ^, centered at 
the origin, is ruled out by present LEP bounds on chargino masses |]^ . As the dependence 
of the rate ( ]25|) on the SU{2)l gaugino mass M is essentially flat (not shown), these LEP 
bounds effectively translate in our case into excluding /x < QQGeV. The various spikes 
in these figures refiect the points where the overall numerical contribution from the form 
factors (vertex plus self-energy) cancels out and even changes sign. The actual point 
where this cancellation occurs depends on the particular choice of the parameters. 

The fact that in all these figures the ratio (^) is smaller for the heaviest scalar Higgs 
(if °) is not due to the smaller phase space factor, but to the smallness of the form factors. 
We can see in Fig. ^ that in fact B{t — >■ cH^) grows with the A^ (and thus the if°) mass, 
until the phase space begins to close. 

We conclude that a typical value of the ratio (|25|), at large 30 < tan/3 < 50 and for a 
SUSY spectrum around 200 GeV ^ reads roughly 

^SUSY-EW(^^^;^)^lQ-8 ^ (27) 



provided Ma < 120 — 130 GeV. This is larger than the previous reported ratios |2^ by 2 
orders of magnitude, specially in the A^ channel, and it is at least 5 orders of magnitude 
larger than the SM rate B{t — > cHsm) B |13- I^ ia.ct, this feature can be explicitly 
checked in our framework in the limit Mj^o ^ oo in which the CP-even Higgs boson 
h'^ behaves like the SM Higgs boson Hsm- This is already seen in part in Fig.^. By 
further sending the SUSY masses to infinity we indeed recover a very poor FCNC rate 
for /i° ~ Hsm which goes down 10"^^ [g, |TD[. 

Turning now to the analysis of the SUSY-QCD effects, it is clear that they hinge to a 
great extent on the values of the flavor mixing coefficients 5ij. The latter are constrained by 
low-energy data on FCNC. The bounds have been computed using some approximations, 
so they must be taken as order of magnitude limits rather than as accurate numbers. 
They read as follows [^ p9|] : 



I (^12 1 < ■l^/m^^nc/bQQGeV , 

\5i2.\ < .mS^/m^mi/bmGeV , 

\52z\ < S.2m~^mi/{hmGeVf . (28) 

In using these bounds we make use of SU(2) gauge invariance to transfer the experimental 
information known from the down-quark sector (for example from BR{b — i> s 7), where 
the bound on ^23 is obtained) to the up-quark sector. It means that after soft SUSY 
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breaking, but before SSB, the LL blocks of the up-squark and down-squark mass matrices 



must satisfy the following relation [^ 



{MI)ll = K {MI)ll K^ , (29) 

where K is the Cabibbo-Kobayashi-Maskawa matrix. Thus if M~ is parametrized as in 



eq. (15), then {M~)ll inherits a similar form with a new set of mixing coefficients 5ij which 



differ from the previous ones by factors of (9(1). 

We use the same input parameters as in the electroweak contributions (p6[) , plus the 
specific parameters of the SUSY-QCD sector, namely the gluino mass rrig and the mixing 



coefficients 5ij, as follows: 



Trig = 200 GeV , 



0.03 0.03 
0.03 0.4 I . (30) 

0.03 0.4 

As for the strong coupling constant we have used the value as{mf) =0.11. A comment is 
in order for the present set of inputs: we have introduced in ([26| ) the lightest stop mass as 
an input, and this stop is mostly a t/j. However, in this new parametrization we introduce 
this mass as the lightest Mq, mass, which again will be mostly a tn. Notice furthermore that 
the chosen entries in the mixing matrix (|30D are moderate in the sense that they do not 
saturate the permissible upper limits ([28|). Here we are using the additional constraint 
that the squark masses (which are obviously affected by the values of the parameters 
Sij) cannot be too light. As mentioned above, although we use a fixed "typical" choice of 
inputs, a systematic scanning has been performed over the parameter space. In particular, 
the maximum values for the rates (to be compared with the typical ones) have also been 
pinned down (see later on). 

In analyzing the SUSY-QCD effects we find that again the largest contribution comes 
from the right-handed form factor in eq. (p^), but this is only because up to now we 
have chosen not to introduce mixing between right-handed squarks. We have plotted 
the evolution of the ratio (|25| ) with some parameters of the MSSM in Fig. ^ As it 
reads off eq. (0) the most relevant parameter in the SUSY-QCD analysis is the mixing 
coefficient between the 2nd and 3rd generation of LH squarks, which is the less restricted 
one of the three in eq. (|2^) . In Fig. ^ it is shown that by changing ^23 by 3 orders 
of magnitude the ratio ( P3| ) increases by 6 orders of magnitude! - a fact that can be 
traced to the quadratic dependence on the mixing coefficient. Worth noticing in Fig. Qd 
is the fi parameter dependence of the SUSY-QCD effects, which enters through the qiQuh 
coupling. The ratio (p5|) can be pushed up to values of 10~^ irrespective of the sign of fi. 
Again the central region of /j, is excluded by present LEP bounds on the chargino masses. 

The evolution with the gluino mass (Cf. Fig. ^c) is asymptotically quite stable, showing 
a slow decoupling. Thus, even for gluinos as heavy as 500 GeV the rate for the top quark 
decay into the lightest CP-even Higgs boson (t -^ ch^) can stay above 10~^. Finally 
in Fig.|^d we have plotted the evolution with the pseudoscalar Higgs mass. It is seen 
a behavior similar to the EW case (Cf. Fig. ^) although scaled up a factor 10^ — 10'^. 
Obviously, the one- loop Higgs mass relations play an essential role here (missed in ||24|| ) 
in that the h^ is bound to have a mass (below ~ 135 GeV) which is higher than in the 
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tree-level case (below Mz)- Still the asymptotic behavior of B{t — ^ ch^) versus M^o is 
sustained a long while around 10~^ for a SUSY spectrum of a few hundred GeV. 

Figs. |a and b display the change of p5| ) with the lightest stop and sbottom masses 
within the allowed region. It is clear that m^^ plays an important role, which is due to the 
fact that by changing this parameter one is also changing the value of the mixing angle 
between the LH and RH stops. On the other hand the actual value of mi is seen not to 
be that important in the A^ and h^ channels. 

In Fig.^c we plot the ratio (p5|) versus tan/3. We see that the evolution of the SUSY- 



QCD contributions as a function of tan/5 is quite stable for the h^ and A^ channels, and 
so it does not matter much, in the SUSY-QCD case, whether we are in the low or in 
the high tan/3 region. The pronounced downwards spike around tan/5 ~ 10 for the if" 
channel is due to the change of sign of the form factor. 

For completeness we have explicitly probed the impact on the FCNC rates in the 
presence of intergenerational mass-mixing terms in the right-handed sector of the model. 
In order to maintain the number of parameters under control we have set {Sij)fiji = {6ij)LL. 
We find that the inclusion of the additional mixing coefficients does increase the FCNC 
rate, but we do not plot the result since it just amounts to a total contribution which 
is at most twice the old result with {5ij)RR = 0. And this feature holds for any set of 
input parameters ( p6| , pOD . Therefore, we conclude that the old coefficients {5ij)LL alone 



-the only ones that are naturally generated within RG-based models ^\ - already bring 
about the bulk of the FCNC rates (|2S|). 

We wish to emphasize that although we have used a common subset of inputs -see 
eq. (^)- to compute the SUSY-EW and the SUSY-QCD loops, the enhancement sources 
for the two types of contributions are entirely different. Thus, whereas the SUSY-EW 
effects are much sensitive to extreme values of the parameter tan/3 through the Yukawa 
couplings (y), the SUSY-QCD effects (which are the leading ones in our calculation) are 
not particularly sensitive to tan/3, as confirmed in Fig. ^. Instead, they are basically 



dependent on the experimentally allowed values of the flavor mixing coefficients (15), as 
it is plain in Fig. ^. 

From the previous numerical analysis we confirm that the preliminary MSSM treat- 
ment of the FCNC decay of the top quark into Higgs bosons |2^ was fairly incomplete 



since important effects from Higgs particles in the loops were not included, and moreover 
the qiQRh vertices were not taken into account. As a result the potentially large contri- 
butions coming from the trilinear soft SUSY-breaking terms At^b, and from the higgsino 
mass parameter fi were missed. Moreover, the pattern of quantum effects in the three 
Higgs channels is affected in an essential way by the one-loop Higgs mass relations of the 
MSSM. From our rigorous computation we have been able to show that B{t -^ ch) in the 
MSSM can typically be of order 10"*^ for the electroweak contributions, and reach 10~^ 
for the QCD contributions. In some of the channels this amounts to having rates that 
are one to two orders of magnitude larger than previous estimates. For an assessment of 
the impact of these results on experimental searches, and a discussion of the maximum 
attainable rates, see Sect. 6. 
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5 The decay t ^ c g 

We have already mentioned that the FCNC decays of the top quark into gauge bosons, 
t ^ c V, may receive important contributions from SUSY physics, and this has been 
confirmed by exphcit calculations in the literature |jT^-|l^. In particular, the situation 
with the decays t ^ cV into electroweak gauge bosons V = j, Z is that, in spite of the 
SUSY enhancements, the final rates are far insufficient to be seen, except for very especial 
circumstances associated to the possibility of wave-function threshold effects [l^. These 



effects, however, are "point-like coincidences", so to speak, which we deem to be very 
unlikely. Indeed, a similar situation occurs e.g. when computing the SUSY corrections 
to the conventional decay of the top quark, t -^ W^h |^. Also in this case there are 
particular combinations of sparticle masses which fall within the (very) narrow range 
where wave-function threshold enhancements occurQ. To be sure, this circumstance should 
be considered exceptional and so the fairest conclusion ought rather to be that the SUSY 
effects on t — * W^h are generally small, and similarly that the SUSY rates for the FCNC 
decays t — >■ c (7, Z) are far below the experimental possibilities. 

In the specific case of the gluon channel t ^ c g one could also argue that it could 
be enhanced by threshold effects, only if the top quark, gluino and stop masses turned 
out to satisfy the peculiar relation nit — ""^7+ f^g- Nonetheless a relation like this is 
not only contrived, but it is already ruled out by the current bounds on gluino and stop 
pl[| : m~ > 180 GeV, m^ > 80 GeV. Fortunately, in contrast to the electroweak 



masses 



gauge boson channels just mentioned, alternative important corrections to t — >■ c (7 are 
possible from gluino-mediated FCNC loops. This subject has recently been addressed 
in Ref. [jl^], but due to discrepancies of this reference with previous calculations by the 



authors of Refs. |I^, |T5| we wish to reanalyze this decay in our framework. In this way 
we hope to further clarify the situation and at the same time to use t — > c (? as a fiducial 
observable with which to better compare the results that we find for the t ^ c h channels 
within one and the same set of assumptions. 

The diagrams contributing to t ^ c g are similar to those in Figs.|I|a,c and ^ after 
replacing the Higgs boson h with the gluon g. The corresponding SM contribution involves 
the W gauge boson and the bottom quark fiowing in the loop. However, in the conditions 
of the foregoing study of t — *> c /i, we have seen that the set of electroweak diagrams in 
Figs. Ill-^ furnish in general a negligible contribution as compared to the gluino-mediated 
contributions in Fig.|^ - and we find that this does not change for the decay t —>■ c g. 
Therefore, we only report on the corresponding SUSY-QCD effects. 

A crucial point in computing this type of effects in the case of t — > c /i was to realize 
that there are diagrams in the electroweak eigenstate-basis with an helicity flip at the 
gluino line. As a result the decay t ^ c h received contributions "proportional" to the 
gluino mass -Cf . eq. ( p3D . These terms cause the corresponding decay rate to fall off very 
slowly with m~, and they even produce a local maximum with respect to this parameter 
(recall Fig.^). Similarly, it turns out that the decay t ^ c g can also have this kind 
of enhancements, although from a slightly different origin (see below). This fact, with 



which we agree with the recent calculation of Ref. |T^ , was missed in Refs. [p!4| , |T^ and it 
led them to speculate on the existence of additional gluino-mediated FCNC interactions 
in the RH sector in order to reach higher rates. Although these additional terms are 
possible, in principle, they are unnatural in a RG-based framework; and what is more. 



^See Ref. g for details, and in particular Figs. 6, 7 and 8 of that reference. 
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they are actually unnecessary to potentially achieve the desired enhancements. 

The analytical calculations ioT t ^ c g are similar to those for t ^ c h, and so we just 
report on the numerical analysis, which we present in Figs. ^y. In all these figures we 
exhibit the evolution of the ratio B{t -^ c g) - defined as in eq. (pS]) with g replacing h - for 
two different scenarios, namely, when only mixing between LH squarks is allowed (solid 
lines), and when mixing mass terms are allowed also in the RH squark sector (dashed 
lines). If one sticks to the more conservative FCNC gluino-mediated interactions in the 
pure LH sector, we see that one may reach important enhancements which can bring 
B{t ^ c g) up to 10^^-10^^ for allowed sparticle massesQ. However, as in the t —>■ c h 
case, we find that the inclusion of FCNC in the RH sector do not modify the order of 
magnitude of the results; in fact the RH interactions only enhance the previous result by 
at most a factor of 2. 

In the calculation, when including the RH terms, we have set {Sij)jiR = {Sij)LL and 
used the fact that the bounds on (%)ll are interwoven with the bounds on the squark 
masses |2^ as dictated by eq. ([T5|), similarly as we did in our analysis of the t ^ c h 
decays in Sect. 3. 

In Fig.|^ we display the evolution of B(t ^ eg) as a function of the mixing parameter 
^23 and the gluino mass nig. We see that Bit — > eg) is also proportional to (^23)^ 
as B(t -^ eh). The evolution with the gluino mass also exhibits a local maximum, 
characteristic of the terms "proportional" to rrig that are triggered by the helicity flip in 
the gluino line. As we said before, this flip occurs even if we only allow intergenerational 
mass-mixing terms in the LL block of the mass matrix - a feature which was completely 
overlooked in Ref. ||15||. The origin of the flip in the electroweak-eigenstate basis is as 
follows. In the one-loop vertex diagram there exists a FCNC gluino interaction with 
the LH charm quark and the LH stop, where the latter stems from a RH stop that 
has mutated (through a mixed propagator) into a LH state. This situation is favorable 
because of large LR mixing in the stop sector. Schematically, one can think in terms of the 
following picture: from the gluon vertex in the loop there emerges a tn that subsequently 
undergoes the series of transitions in ^ ii —* ei |. As the structure of this loop enforces 
a mass insertion in the gluino propagator, it leads to an enhancement of B{t —* eg) which 
nevertheless falls off with the gluino mass much faster than in the B{t ^ eh) case (Cf. 
Fig.|c). 

The evolution with the electroweak parameters tan/? and fi is presented in Fig.|^. In 
this process they only enter as inputs in the squark mass matrix, but not in the couplings, 
so the dependence of B{t ^ eg) on them is rather mild. In the presence of mixing in both 
the LH and RH sectors, the physical squarks masses are smaller than the corresponding 
ones with mixing only in the LH sector. Eventually these particles can be lighter than 
present bounds, thus the parameter space is further constrained. This is shown in Fig.^] 
where the dashed lines (corresponding to 6ll = ^rr) are cut off at points in which the 
present bounds on squarks masses would be violated. Finally, we report in Fig.|^ on the 
evolution with the lightest stop mass and the soft SUSY-breaking trilinear coupling Af. 
The shaded region in Fig. ^d is excluded as it incompatible with our choice of parameters, 

^When comparing with Ref. |g^, in the unconstrained MSSM case, we point out that our rates are 
smaller only because we use an updated (larger) value of the gluino mass and also a larger value for the 
lightest stop mass - see eq. (|2^)- which is the most sensitive one. 

^The difference with the t ^ ch case lies in the fact that the first transition in -^ II is already 
possible at the Higgs vertex II in h in the loop, but it cannot take place at the gluon vertex. 
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and the various lines end up where the large mixing would induce too light squark masses. 
The rate is seen to be sensitive to both parameters m^^ and At, but the actually permitted 
region for the latter is quite narrow. 

To summarize, upon comparing the decay t ^ c g with the relevant decays t ^ c h 
under study we find that the latter can be relatively very important, and most likely they 
are the dominant FCNC top quark decays in the MSSM. As a matter of fact, whereas the 
critical level 10~^ is never surpassed by the gluon channel - even with LH+RH effects and 
for sparticle masses within the current limits-, the t —>■ c h^ and t —>■ c A^ modes, instead, 
may well crossover the visible level already with the more conventional LH contributions. 
The maximum values for all these decays are discussed in the next section, where we also 
present our conclusions. 

6 Discussion and conclusions 

We have studied the SUSY-EW and SUSY-QCD contributions to the leading FCNC top 
quark decays. We have mainly concentrated on the top quark FCNC decays into the Higgs 
bosons of the MSSM, namely t ^ ch {h = h^, H^, A^), using a mass spectrum motivated, 
but not fully restricted, by model building and Renormalization Group Equations. And 
we have found that a full treatment of the SUSY-QCD contributions may greatly enhance 
the FCNC width by some orders of magnitude with respect to the FCNC decay rates into 
gauge bosons t -^ cV {V = 'y,Z,g). Due to the crucial role played by the SUSY-QCD 
effects, we have reconsidered the FCNC decay into gluons, t ^ eg, which is the most 
promising one among the FCNC top quark decays into gauge bosons. And we have shown 
that under similar assumptions of gluino-mediated FCNC couplings in the LH sector, it 
could be enhanced up to near the visible level, but below the rates of the Higgs channels. 
The additional RH interactions, if present at all, could produce a further increase of the 
rates, but it just amounts a factor of 2 at most. The remaining gauge boson channels 
give maximum rates in the MSSM that are far below experimental possibilities, except in 
highly unlikely circumstances. 

With a SUSY mass spectrum around 200 GeV, which is above the current absolute 



LEP bounds ||3T|, the different contributions to the Higgs channels are typically of the 
order 

^SUSY~QCD(^^^^) ^ ^Q-5_ ^3^) 

However, by stretching out a bit more the range of parameters one can reach (for some 
of the decays) 

^susY-QCD^^^^^) ^ 5x10-^ (32) 

The difference of at least two orders of magnitude between the SUSY-EW and SUSY- 
QCD contributions makes unnecessary to compute the interference terms between the two 
sets of amphtudes, but if the limits on ^23 (eq. (|28|)) become eventually more restrictive 
then they should be taken into account. 

We have obtained the maximum rates (E^ from a general search in the MSSM pa- 



rameter space within the 1 TeV mass region. In Fig. |10| we present the maximum values 
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that can be reached hj B{t —>■ cX) for each of the processes presented. In Figs.lTOa and b 



we show the maximized B(t ^ ch) as a. function of the pseudoscalar Higgs boson mass 
by taking into account only the SUSY-EW contributions and the SUSY-QCD contribu- 
tions respectively. We have performed a systematic scanning of the parameter space of 
the MSSM, with the various masses constrained between the present exclusion bounds 
and the ITeV upper bound, and using the other bounds from Sect.||- in particular, see 
eq. (^). Perhaps the most noticeable result is that the decay into the lightest MSSM 
Higgs boson (t —>■ ch^) is the one that can be maximally enhanced and reaching values 
of order B{t —>■ ch^) ~ 10^'^ that stay fairly stable all over the parameter space. The 
reason for this dominance is that the decay t ^ c/i" is the one which is more sensitive 
to At, a parameter whose natural range reaches up to about ITeV. For more moderate 
values of At, however, the dominant Higgs decay mode can be t ^ cH^, but in this case 
the corresponding rate (also that of t — ;> cA'^) undergoes a rapid fall-off with M40 (Cf. 
Fig. pUjb) as it was also the case in the unmaximized situation in Fig.^. In general, the 
(next-to-leading) dominance oi t ^ cH^ in Figs.p!0|a and b is confined to a small corner 



of the parameter space which we pick up during the process of maximization. Quite in 
contrast, the FCNC top quark decay into the lightest Higgs scalar can have an observable 
ratio in a large portion of the parameter space, and in particular for almost all the range 
of Higgs boson masses. Needless to say, not all of the maxima can be simultaneously 
attained as they are obtained for different values of the parameters. 

Moreover, we have also found (Cf. Fig.pHjc) the maximum FCNC rate of the gluon 
channel in the MSSM. Under the RG-based assumption of only mixing in the LH sector, 

5SUSY-QCD(^_^^^)< ^Q-5^ (33) 

but it never really reaches the critical value 10~^, which can be considered as the visible 
threshold for the next generation of colliders (see below). The visible limit for the gluon 
channel can only be picked up if one maximizes the ratio under the assumption of both 
LH and RH similar contributions, but even in this case the limit is only barely reached. 
However, we emphasize that the right order of magnitude can already be achieved with 
only flavor mixing in the LH sector. So for this decay we confirm the recent analysis 
of Ref.|]2B| in contrast to that of Ref. [^. In Fig.pHjc we have plotted the maximum 
value oi B(t -^ eg) as a. function of 623 after scanning for the rest of the MSSM mass 
parameters within the 1 TeV range. Interestingly enough we see that the maximum rate 
is not reached for the highest possible value of the flavor mixing parameter 623, the reason 
being that the physical squark masses - obtained after diagonalizing the 6x6 squark 
mass matrix in (flavor) x (chiral) space - are constrained from below, so that the higher is 
623 the smaller is the maximally allowed value for the left-right mixing in the stop sector 
-a result which also holds for the case of the Higgs channels (p . 

The remarkable enhancements obtained for the top quark FCNC Higgs decays (0) 
could bring them to detectable levels in large portions of the parameter space and not 
just in small optimized domains. To assess the discovery reach of the FCNC top quark 
decays in the next generation of accelerators we take as a guide the estimations that have 
been made for gauge boson final estates ||3^. Using the information mentioned in Sect.|l| 



and assuming that all the FCNC decays t — > c X {X = V, h) can be treated similarly, we 
roughly estimate the following sensitivities for 100 fb~^ of integrated luminosity: 

LHC :S(t^cX) > 5x10"^ 
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LC:5(t^cX) > 5x10"^ 
TEV33 -.Bit ^cX) > 5 X 10"^ . (34) 

Therefore, LHC seems to be the most suitable colhder where to test this kind of phenom- 
ena. The LC is hmited by statistics (due to much smaller top quark cross-section) but in 
compensation every collected event is clear-cut. So this machine could eventually be of 
much help, especially if we take into account that it could deliver 500 fb~^ per year 
during the highest luminosity phase. In fact, a high luminosity e~^e~ super-collider offers 
the greatest potential for high precision top quark physics and it constitutes an ideal com- 
plement to concomitant LHC experiments. So, even if no SUSY particle is seen at LEP II, 
there is indeed a possibility to pin down quite a few FCNC t ^ ch and/or t ^ eg decays 
at the LHC first, and later on to perform a more detailed study at the LC. The situation 
with the Tevatron, as we see, is much more gloomy, since the required FCNC rates of 
~ 10^^ cannot be attained unless we artificially search in some remote (unnatural) corner 
of the parameter space. Hence, the Tevatron (even TeV33) seems to be of no much help 
in this regard, unfortunately. 

Of course, with the range of SUSY masses that we have used in our analysis, one 
could try to detect some of the sparticles directly at the super-colliders. But this is al- 
ways harder than dealing with conventional particles, and in any case the FCNC method 
is complementary to direct searches. Therefore, from the enormous experience gained in 
the study of top quarks, a dedicated analysis of the potential few hundred FCNC top 
quark decays per year at the super-collider LHC could provide a more clear identification 
of physical effects beyond the SM. This is so not only because the SUSY-FCNC rates could 
well be within the experimentally feasible ranges (|5^), as we have seen in our numerical 
analysis, but also because the signatures are likely to be separable from the background. 
For instance, at high tan/? > 30 all of the final state Higgs bosons h = /i°, H^, A° would 
mainly disintegrate into bb pairs, h —>■ bb. These decays are also highly augmented as com- 
pared to the SM prediction, and include large MSSM radiative corrections |^4[. Although 
there is an interval of Higgs masses where the alternative decay h -^ ZZ is possible and 
sizeable for an SM Higgs, this decay mode is not dominant for the SUSY Higgs bosons 
h = H^, A^ in the relevant region, and it is never kinematically accessible to h^. In this 
region we expect that the FCNC final states ti ^ ch, ch are to be observed mainly as 
ebb, ebb, and they should be effectively tagged through high pt charm-quark jets and 
large invariant mass for the recoiling bb pairs. On the other hand, for moderately small 
tan/? > 1 the situation is more cumbersome since the gluino-mediated FCNC rates for 
t ^ eh can be equally sizeable as for high tan/3, but in turn the decay pattern of the 
neutral Higgs bosons is much more complicated, and this is so the heavier are the Higgs 



bosons 1^. Notwithstanding, since in our case the Higgs bosons must satisfy nih < rrit, 
it turns out that both the lightest CP-even and the CP-odd states h = h^,A'^ still do 
preferentially decay into bb pairs. But not necessarily so the heavy CP-even Higgs H^, for 
in this region it first cascades down predominantly (90%) into h^ h^ (except in a narrow 
interval centered around Mj^o = 130 GeV where if ° — > bb again dominates because the 
H^h^h^ -coupling changes sign); and as a result the observable bb pairs in the final state 
(emerging from the two secondary - real or virtual - decays h^ -^ bb) are much softer 
than in the previous situation. Fortunately, we have seen that among the FCNC decays 
t ^ eh the one corresponding to h = H^ is typically the most suppressed one, except 
in some corners of parameter space. In contrast, the most relevant decay (both under 



typical and optimized conditions) is that oi h = h^. The corresponding channel is always 
kinematically open in the MSSM and it can be tagged through the final state signature 
ebb mentioned above, which is dominant for any tan/3 > 1. As for the FCNC decays in 
the gluon channel, t ^ eg, here also a high pt charm-quark jet against a highly energetic 
(~ mj/2) gluon jet is the indelible imprint of this decay. Altogether these features should 
be very helpful to mostly reduce all backgrounds and provide a clear-cut signal of rare 
top quark decays beyond the SM. 

To conclude, the FCNC decays of the top quark are such rare events in the SM 
(especially the top quark decay into the Higgs boson) that their observation at detectable 
levels should be interpreted as an extremely robust indication of new physics. At present 
the most popular (phenomenologically compatible ^3\ and technically fully consistent 0) 



quantum field theoretical extension of the SM is the MSSM; therefore, the effective tagging 
of FCNC top quark decays in the major accelerators round the corner could well be a first 
step into discovering SUSY. From our analysis we find that around the loci of maxima of 
the rates for the leading FCNC top quark decays, t ^ eh and t ^ eg, we roughly have 
the following situation (under the assumption of fiavor mixing only in the LH sector): 

5 X 10-6 < B{t -^ eg)^,^ < B{t -^ eh)^.. < 5 x 10'^ (35) 

In both types of decays the dominant effects come from SUSY-QCD. However, it should 
not be undervalued the fact that the maximum electroweak rates ioi t ^ eh can reach 
the 10~6 level, which amounts not only to saying that they are 7 orders of magnitude 
greater than the maximum SM ones, but also to noticing that they are on the verge of 
being detectable. Last but not least, we stress once again that the largest FCNC rate 
both from SUSY-QCD and SUSY-EW is precisely that of the lightest CP-even state, a 
very fortunate fact which can hardly be overemphasized as t -^ c/i° is the only Higgs 
channel that is phase-space available across the whole MSSM parameter space. 

At the end of the day, we should seriously bear in mind the possibility of observing 
FCNC decays of the top quark, and so we ought to be prepared to collect a few hundred, 
perhaps a few thousand, events of this sort at the LHC, together with a few dozen (al- 
though crystal- clear) t ^ eX {X = h,g) decays at the LC. If that would be the case, 
we could eventually find ourselves in the process of discovering SUSY dynamics at work 
before being able to directly produce (or clearly identify) any supersymmetric particle. 
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Figure Captions 

Fig. 1 One-loop SUSY-EW vertex diagrams for tlie decay t -^ ch {h = h^,H^,A^). Here 
d {d{a,b}) represent mass-eigenstate down type quarks (squarks) of any generation. 

Fig. 2 One-loop SUSY-EW diagrams contributing to the mixed t — c self-energy, with a 
notation similar to that of Fig. |l|. 

Fig. 3 One-loop SUSY-QCD diagrams for the decay t ^ ch: (a) vertex diagram, (b) 
mixed t — c self-energy. tt{a,/3} stand for mass-eigenstate up-type squarks of any 
generation. 



Fig. 4 Evolution of the SUSY-EW contributions to the ratio (|25|) with (a) tan/3, (b) the 
trilinear coupling A^, (c) the higgsino mass parameter /x, and (d) the pseudoscalar 
Higgs mass M^o. The rest of inputs are given in eq. (p6|). 



Fig. 5 Evolution of the SUSY-QCD contributions to the ratio (^) with (a) the mixing 
parameter ^23 between the 2nd and 3rd squark generations, (b) the higgsino mass 
parameter /x, (c) the gluino mass nig, and (d) the pseudoscalar Higgs mass M^o. 
The rest of inputs are given in eqs. (PB|) and 



Fig. 6 Evolution of the SUSY-QCD contributions to the ratio ( ^5]) with (a) the lightest 
sbottom mass (m^J, (b) the lightest stop mass (mj-J, and (c) tan/3. 

Fig. 7 Evolution of the SUSY-QCD effects on B(t —>■ eg) as a function of (a) the mixing 
parameter ^23, and (b) the gluino mass rrig. The results are shown under the 
assumptions of: mixing only in the left-handed squark sector (solid); and equal 
mixing in the left- and right-handed squark sectors (dashed). 

Fig. 8 As in Fig.|^ but as a function of (a) the higgsino mass parameter fi, and (b) tan/3. 

Fig. 9 As in Fig.|^ but as a function of (a) the lightest stop mass rra^^ and (b) the soft 
SUSY-breaking trilinear top-squark coupling At. 

Fig. 10 (a) Maximum value of B{t —>■ ch), obtained by taking into account only the SUSY- 
EW contributions, as a function of M^o ; (b) as in (a) but taking into account only 
the SUSY-QCD contributions; and (c) maximum value of B(t ^ eg) as a function 
of the intergenerational mixing parameter ^23 in the LH sector. In all cases the 
scanning for the rest of parameters of the MSSM has been performed within the 
phenomenologically allowed region. 
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